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^ . ABSTRACT 

I ^ I ' We present the results of numerical simulations of continuum-driven winds of 

stars that exceed the Eddington limit and compare these against predictions 
, from earlier analytical solutions. Our models are based on the assumption 

CO . that the stellar atmosphere consists of clumped matter, where the Individ- 

' ual clumps have a much larger optical thickness than the matter between the 

clumps. This 'porosity' of the stellar atmosphere reduces the coupling between 
OO , radiation and matter, since photons tend to escape through the more tenuous 



gas between the clumps. This allows a star that formally exceeds the Edding- 
ton limit to remain stable, yet produce a steady outflow from the region where 



X 

5h ' the clumps become optically thin. We have made a parameter study of wind 



models for a variety of input conditions in order to explore the properties of 
continuum- driven winds. 

The results show that the numerical simulations reproduce quite closely 
the analytical scalings. The mass loss rates produced in our models are much 
larger than can be achieved by line driving. This makes continuum driving 
a good mechanism to explain the large mass loss and flow speeds of giant 
outbursts, as observed in 77 Carinae and other luminous blue variable (LBV) 
stars. Continuum driving may also be important in population III stars, since 
line driving becomes ineffective at low metalicities. We also explore the effect 
of photon tiring and the limits it places on the wind parameters. 
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1 INTRODUCTION 

Massive, hot stars continuously lose mass through radiation driving. The most commonly 
explored mechanism is line driving, wherein the scattering of photons by ions in the stel- 
lar atmosphere transfers momentum from the radiation field to the gas. This mechanism 
results in a quiescent m ass loss with mass loss rates ranging up to about 10^^ Mq yr~^ 
( ISmith fc Owockill2006h . which can explain the winds of most r nassive, hot stars. In fa ct, for 
most stars, the observed mass loss rates are considerably less ( IVink fc de Koteiil2005l ). 

However, some stars, most notably luminous blue variable stars (LBVs) such as r] Cari- 
nae, experience outbursts with mass loss rates several orders of magnitude higher than can be 



expla ined through line driving (jPavidson fc Humphreys 



1997; Smith 2002: Owocki fc van Mark 



20081 ). In the case of r] Carinae, the 1840's outburst is inferred to res ulted in the election of ca. 



10-20 Mq over a time lasting several years, even up to a decade (jPavidson &: Humphreys 



19971 ). While short compared to evolutionary timescale of millions of years, this is much 
longer than the typical dynamical timescale of hours, characterized by either free-fall time 
or interior sound travel time across a stellar radius. Thus in contrast to SN "explosions" that 
are effectively driven by the overpressure of superheated gas in the deep interior, explaining 
LBV outbursts requires a more sustained mechanism that can drive a quasi-steady mass 
loss, a stellar wind, from near the stellar surface. In the case of rj Carinae the outburst was 
accompanied by a strong increase in radiative luminosity, very likely making it well above 
the Eddington limit for which continuum driving by just electron scattering would exceed 
the stellar gravity. The reason for this extended increase in luminosity is not yet understood, 
and likely involves interior processes beyond the scope of this paper. Instead, the focus here 
is on the way such continuum driving can result in a sustained mass loss that greatly exceeds 
what is possible through line opacity. 

A key feature of continuum driving is that, unlike line driving, it does not become 
saturated from self-absorption effects in a very dense, optically thick region. Indeed, since 
both continuum acceleration and gravity scale with the inverse square of the radius, a star 
that exceeds the Eddington limit formally becomes gravitationally unbound not only at the 
surface but throughout. Clearly, this is in contradiction to the steady surface wind mass 
loss observed for these stars. N.B. Contrary to what is sometimes claimed, this would not 
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automatically destroy the entire star. Although radiative acceleretion might overcome gravity 
locally, the total energy in the radiation field would not suffice to drive the entire envelope 
of the star to infinity (this is known as 'photon tiring'). Instead, the outward motion of the 
gas would quickly stagnate and matter would start to fall back. Nevertheless, the net result 
would not resemble a steady wind. 

This problem can be resolved by assuming that the stellar material is clumped rather than 
homogeneous, with the individual clumps being optically thick - and therefore self-shielding 
from the radiation - whereas the medium in between the clumps is relatively tranparent 
to radiation. T his so called ' poros ity effect' can lead to a reduced coupling between matter 



and radiation (Shaviv 



1998 



2OOOI ). The photons tend to escape through the optically thin 



material between the clumps without interacting with the matter inside the clumps. This 
implies that a star that formally exceeds the Eddington limit can remain gravitationally 
bound and would only exceed the effective Eddington limit at the radius where the individual 

clumps themselves become optically thin. 

The structure of such a star should therefore look as foll ows Jshavivlknnibi. Deep inside 



the super-Eddington star, convection is necessarily excited (jjoss. Salpeter fc Ostrikeiill973l ) 
such that the radiation field remains sub-Eddington through most of the stellar interior. 
At low enough densities where maximally efficient convection cannot sufficiently reduce 



the radiative fiux, the near 



possible instabilities (lArons 



ddington luminosit y necessarily excites at least one of several 



1992 



Shaviv 



2001al . e.g.,) which give rise to a reduced opacity. 



This 'porous' layer has a reduced effective opacity and an increased effective Eddington 
luminosity. Thus, the layer remains gravitationally bound to the star. At lower densities still, 
the dense clumps become optically thin and the effective opacity approaches its microscopic 
value. From this radius outwards, the matter is gravitationally unbound, and is part of a 
continuum-driven wind. 

A detailed a nalytical study of this paradigm was carried out by Owocki, Gayley & Shaviv 
(2004), hereafter IOGSO4J . This predicted that continuum-driven winds can produce high mass 
loss rates (^ 10~^ Mq yr~^) at intermediate wind velocities (10^ — 10'^ km s~^). Here we 
test these analytical predictions with numerical simulations of winds from super-Eddington 
stars. 

In addition to LBVs, which are the specific objects we study here, other types of astro- 
nomical objects can exceed the Eddington limit and the refore experien ce similar continuum- 
driven winds. These include for example classical novae (jShavivll2001bl ) or high accretion rate 
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Figure 1. Left: The effect of continuum driving on a sub-Eddington, line-driven wind for Fe = 0.1. The solid lines depict 
the wind parameters for the CAK force only, the dashed lines show the joint effect of both CAK line driving and continuum 
scattering. Clearly, the influence of the continuum term is negligeable. Right: The same variables , but for Fe = 0.5. Here the 
difference caused by continuum driving is significantly more pronounced. 



accretion disks around black holes (IBegelmaru 120061 ) . In fact, classical nova eruptions clearly 



exhibit steady continuum-driven winds, indicating that the mass loss rate is somehow being 
regulated and likely to be described by the porosity model and the same continuum-driven 
wind analyzed here. 

The layout of this paper is as follows. In ||2] we show the effect of continuum scattering 
on a su b-Eddington, line-driven wind. In fj3]we summarize the analytic results obtained by 
OGS04. In ® we describe the numerical methods that we have used for our simulations. 



^ shows the results of our simulations and the comparison with the analytical predictions. 
In f|6] we discuss the effect of photon tiring and show how it influences the results of our 
calculations. Finally, in f|7] we end with a summary and a discussion. 



2 SUB-EDDINGTON LIMIT 



As long as a star does not exceed the Eddington limit, continuum interaction alone can- 
not drive a stellar wind. However, it does influence the parameters of a (primarily) line- 
driven wind. Figure [1] demonstrates the influence of continuum driving on the stellar wind 
parameters. For the s e sira ulations we use the CAK line driving formalism described in 
ud-Doula fc Owockil (120021 ) . including both the finite dis k correction factor and the modi- 
fied gravity scaling of the line force (eq. [19] of |OGS04| ). For the same stellar parameters 
(50 Mq, 25 Rq, 50 000 K) and CAK parameters (Q = 2.0 x 10^ a = 2/3) but two different 
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Eddington parameters ( Fg = L^^/L-^dd = 0.1 and 0.5,) we simulate the stellar wind with and 
without the continuum driving term. 

Note that unlike the purely continuum-driven winds in the next section, these simulations 
do not contain any porosity effect. Namely, the continuum driving term in the equation of 
motion is simply: 

S'cont — '^e ~j 

with the continuum opacity set to 0.3 g/cm^. These parameters give us an Eddington 
luminosity of -Lsdd = 4:TtGcM/ = 8.33 x 10^^ erg/s. Figure [1] shows both the density 
(relative to the density at the stellar surface) and wind velocity as a function of radius for 
these two simulations. Clearly, for a low Eddington factor, the effect of continuum interaction 
on the wind parameters is negligeable. However, once the star approaches the Eddington 
limit, continuum driving increases the mass loss rate and decreases the terminal velocity. 
For Fe = 0.1 the mass loss rate increases from 1.80 x 1O~^M0 yr~^ to 1.92 x 1O~^M0 yr~^ 
because of continuum driving. For Fg = 0.5, these numbers are 2.18 x 1O~^M0 yr~^ and 
3.22 X IQ-^ Mq yr~^ respectively. 



3 ANALYTICAL APPROXIMATION 



3.1 Wind velocity 

For a radiatively driven outflow from a stellar surface, irrespective of the specific driving 
mechanism involved, the equation of motion can be written as 

(2) 



1-^ 



dv 



dr 



GM^ 2a^ da^ 

+ fl'rad H 3—, 

r dr 



with a the isothermal sound speed, (/rad the radiative acceleration and GM^^/r"^ the inward 
gravitational acceleration. The last two terms on the right hand side are the result of the 
gas pressure gradient. In most cases they be neglected as their contribution to th e velocit y 
is usually negligeable next to the gravitational and radiative acceleration terms ( OGSO4I ). 
This assumption becomes incorrect only in those situations where the terminal velocity of 
the wind is close to the sound speed. If we neglect these gas pressure terms, we can rewrite 
equation ([2]) in a dimensionless form. 



w 



w 



1 + F, 



(3) 



with w = iy/ fesc)^5 Ws = (a/ fesc)^ and Fg = (^rad'^^/GM^. Note that Vesc is defined as 
the escape velocity at the surface of the star i?^ rather than as a local escape velocity. The 
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gravitationally scaled inertial acceleration, 

w' = 



r'^vdv/dr 



(4) 



can be written in terms of the inverse radius coordinate x = 1 — Ri,/r, so that w' = dw/dx. 
This inverse radius coordinate makes for a more practical coordinate system than the radius 
itself, since it scales with the gravitational potential. 

Since the escape velocity of a massive star is typically at least an order of magnitude 
larger than the local isothermal sound sp eed, we can in many cases neglect the bracketed 



term. For the porous atmosphere. 



OGSOJ introduced a dimensionless parameter = i^es/i^, 



which gives the effective opacity of the loc al clump ed medium vs. the opacity of the medium 



if it were completely homogeneous. Like 



OGS04I . we assume that the distribution (/) of 



optical depths of the individual clumps can be described by a truncated power law. 



dr 



Tiar 



-t/tq 



(5) 



with Op ^ being the power-index and tq the optical depth of the thickest clump. Here 

r(ap) is the gamma function, rather than the Eddington parameter. This gives us: 

(l + ro)^---l 



kh 



For the full derivation see 



0GSn4 



(1 - ap)To 



(6) 



As long as the clump characteristics are fixed, the optical depth tq scales with the density 

as, 

^o = p/Po, (7) 
where po is the critical density at which the thickest blob has unit optical depth. This is 
calculated by 

(8) 



Po 



with K, the continuum opacity and ho the porosity length of the thickest clump; the porosity 
length h is defined as L^/l"^, with L being the typical separation between clumps and / the 
size of the clump. 

We thus rewrite eq. ([3]) by substituting, 

Fe FeA;;,(ro), (9) 

for the Eddington parameter Fe = L^/L-^dd- The scaled velocity gradient thus takes the 
form, 

Teh[To{x)] - 1 



w'{x) 



Ws/W 



(10) 
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This is a first-order differential equation that can be integrated using standard numerical 
techniques, starting from the sonic-point initial condition w{x = 0) = Ws, where the vanish- 
ing of both the numerator and denominator on the right-hand side requires application of 
L'Hopital's rule to evaluate the initial gradient, 

w'^ = w'{x = 0) = ^ 



1 + 32ws/ap 



(11) 



By introducing the porosity term into the equation of motion, we have made the ra- 
diative acceleration dependent on the density and therefore on the radius. Under these 
circumstances, the radiative acceleration is only able to drive matter away from the star in 
the outermost layers of the star, producing a steady wind that can last as long as the star 
remains about the Eddington limit. In effect, porosity regulates the mass loss to a level that 
can be sustained in a nearly steady way throughout an LBV outburst. Here again the word 
"outburst" has to be interpreted in context. LBV outbursts can last years, which, though 
short compared to stellar evolutionary time, is quite long compared to a SN explosion. Most 
importantly, this is much longer than a typically djTiamical flow time (ca. a free fall time, 
R/ Vesc, or a wind expansion time, R/voo)- As such, the mass loss can be well modeled in 
terms of steady-state wind solutions that assume the luminosity, etc. are constant over such 
dynamical timescales. 

3.2 Mass loss rate 

The mass loss rate induced by continuum driving is directly related to the luminosity of the 
star and the power index ap. The generic mass loss rate that follows from the porosity is 

M = (12) 
\poJ \hoJ ac 

with H the gravitational scale height of the gas at the sonic point. The sonic point density, 
p*, can be found through condition 

Teh = 1. (13) 

For the simple case ap = 1/2, this has the explicit solution 

-) =4re(re-l) ; ap = l (14) 



which leads to an explicit expression for the mass loss rate (cf. eqn. (77) in |OGS0J) 

1 

with r/o = Iiq/H . For more general cases, |OGS04| also give explicit approximations for the 



M = 4(re-1):-^ ; ap=-, (15) 



expected mass loss rate, but in quoting values for this "analytic porosity model" below, we 
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choose the more accurate approach of solving eqn. (fT3l) imphcitly, and using the resulting 
Pi, to compute the associated mass loss rate from eq. ( |T2|) . 



4 NUMERICAL METHOD 



For our numeric al simulations we use the ZEUS hydrodynamics code (jStone fc Norman 



1992 



Clark 



19961 ). We do our computations on a ID spherical grid with an inflow boundary 



in the center and outflow at the outer boundary. Continuum driving is modeled through an 
accel eration term Vpkh qgrav a dded to the radial component of the equation of motion (see 



also, 



van Marie et al. 



(|2008aD). 

We start our simulation by initializing a ID radial grid with matter moving away from 
the origin with a power law velocity distribution. At the inner boundary, we specify an inflow 
density, higher than the critical density ps to ensure that the sonic point (where the wind 
velocity equals the isothermal sound speed) is inside the grid. This inflow density remains 
constant during the simulation. The inflow velocity is time-dependent, and recalculated 
at each timestep such that the velocity gradient over the boundary is equal to the veloc- 
ity gradient immediately above the boundary radius. Such "floating" boundary conditions 
allow approach to a stable, steady, fl ow solution, and have been succ essfully used for line- 
driven-wind simulations in both on e (lOwocki. Castor fc Rybickilll988l ) and two dimensions 
(jOwocki. Cranmer fc Blondinlll994j ). 

The radial grid is not evenly spaced, but rather chosen in such a fashion that the indi- 
vidual grid size at the sonic point is always smaller than the local scale height. 

For our simulations we assume that the gas is isothermal with a temperature of 50 000 K, 
which implies an isothermal sound speed of approximately 20 km s~^. The opacity of the 
gas is set to 0.4 gr/cm^ (Thompson scattering opacity for pure hydrogen) and r^o fo 1. The 
radius and mass of the central star are set to 50 R© and 50 respectively. A star such 
as this has an Eddington luminosity of 6.24 x 10'^^ erg/s. 

The grid extends from the stellar surface to a distance of ten times the stellar radius. 



5 RESULTS 

A grid of 30 models was calculated, encompassing a parameter space of 0.1 ^ ctp ^ 0.8 and 
1.1 ^ Fg ^ 6.0. Stellar mass, radius surface temperature, opacity and clumping parameter 
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Figure 2. Wind velocity and density as a function of radius for a continuum-driven wind. On the horizontal axis, the local 
radius divided by the sonic point radius. Independent parameters: Ve = 3, Op = 0.5, «)* = 0.001 



?7o were kept constant. A typical result for these calculations is shown in figure [2], which 
shows the density and wind velocity as a function of distance from the star. 
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Figure 3. Wind velocity in units of the sonic point escape radius plotted versus the dimensionless radial coordinate x = 
1 — Rsoniaji'i Comparing results for the analytical solution (dashed curves) and numerical simulation (solid curves) for a 
continuum-driven wind. The upper curves show results for the same parameters as figure [2] with Fe = 3 and a-p = 0.5, while 
the lower curves are for the marginal super-Eddington case Fe = 1.1, with Qp = 0.1. The agreement between numerical and 
analytical solutions is quite good, with only a slightly stronger differences for the Fe = 1.1 case. 
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Table 1. For parameters Fe and given in the first two columns, comparison of the mass loss rates predicted by the analytic 
porosity model (column 3) vs. those obtained in full numerical simulations (column 4), with the fifth column giving the ratio, 
analytical/numerical. The agreement is remarkably good, but with modest, ca. 10% differences for lower ap or F^ near unity. 
The last column gives the ratio between the terminal velocity and the escape speed from the sonic radius for the numerical 
models. 



Table [T] compares mass loss rates for the numerical and analytical solutions for a sample 
of model parameters and Up. As noted above, the analytical mass loss rate is computed 
from eqn. ( fT2i) using the sonic point density derived from implicit solution of eqn. ( |T3l) . 
The comparison shows that the analytical and numerical results coincide very well, with 
maximum differences of about 10% for Fg near unity and small Up. 

The wind velocity can be approximated semi-analytically by integrating equation (fTOl) . 
The results again match the numerical simulation quite well, as plotted in figure [31 

The match is partcularly good for the standard parameter set used in figure Ej with 
Fg = 3 and ap = 0.5. For the marginal super-Eddington model Fg = 1.1 with small power 
index ap = 0.1, which represents an almost pathogical case, the differences are greater, but 
still quite modest. 
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Figure 4. Left: Total mass loss rate in solar masses per year as a function of Fe and ctp. As predicted in the analytical 
approximations (see Q the mass loss rate increases with Fe but decreases with Op. Right: Terminal wind velocity divided by 
escape velocity as a function of Fe and ap. The wind velocity increases with both Fe and Op. The white circles indicate the 
actual datapoints, taken from table [T] 



The complete results of the parameter study are shown in fig. HJ which shows the mass loss 
rate and terminal wind velocity for all simulations. As expected, the mass loss rate increases 
with Fg and decreases with ap as predicted in § [3l The terminal wind velocity increases 
both with Fg and with ap. This is hardly surprising. A high Fg means a strong radiative 
acceleration, which will lead to a higher velocity. A high ap means a weaker coupling between 
matter and radiation. This decreases the mass loss rate, which for a given total available 
energy means an increase in velocity. 

Mass loss rates in general are very high, though one should remember that photon-tiring 
was not included in these simulations. Although photon tiring does not change the mass loss 
rate at the stellar surface, it can prevent that part of the mass escapes the stellar gravity. 
Wind velocities are usually on the order of the escape velocity. 



6 PHOTON TIRING 

One of the most interesting aspects of continuum-driven winds is their hi gh mass loss rates, 



which gives rise to the phenomenon of "photon tired winds" suggested by 



Owocki &: Gayley 



( ]l997l ). That is, winds in which the energy required to drive the gas out of the gravitational 
well is comparable or even larger than the energy available in the radiation field. In §6.11 we 
briefiy cover the analytical description of these winds, and then continue with a numerical 
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study of the case in which there is sufficient energy in the radiation field. The opposite hmit, 
which gives rise to non-stationary solutions will be dealt in the subsequent publication. 



6.1 Analytical approximation 

If the stellar wind is driven by radiation, conservation of energy limits the maximum me- 
chanical luminosity of the star to the radiative luminosity. This effect is known as photon 
tiring and was not included in the simulations described in the previous section. 

The maximum amount of mass that can be driven of by radiation is the mass loss rate 
for which the mechanical luminosity of the wind matches the radiative luminosity, 

^Mv,J ^ Lur = L,, (16) 

with the radiative luminosity of the star. This sets the upper limit for the mass loss rate 
to, 

Mtir = (17) 

"^esc 

Note that this is the same for both line driving and continuum driving, and is unrelated to 
the composition of the stellar atmosphere. It is a fundamental property of the star. From 
this limit we calculate the photon tiring number m, 

M , , 

m = -. — . (18) 

As long as m <^ 1, photon tiring has no significant infiuence on the wind properties, but 
when m < 1 the situation changes. The radiation field of the star becomes depleted as much 
of its energy is used to drive the matter, leaving less energy available to drive the outer 
layers of the wind. The mass loss rate at the surface remains the same since it is set by the 
local radiative acceleration, but the wind velocity will decrease. 

This will also infiuence observations, since radiative energy that is used to drive the wind 
will no longer be part of the visible radiation output of the star. For stars with line-driven 
winds, this is hardly worth considering. Even a Wolf-Rayet star with a powerful wind loses 
only a few percent of its luminosity to the wind. For continuum-driven winds, the effect can 
become quite significant. 

Should m actually exceed unity, the radiation field loses all its energy. As a result, the 
matter that leaves the stellar surface will not be able to reach the escape velocity and start 
to fall back onto the star. This will make it impossible to obt ain a ste ady state solution. 

A analytic analysis of the photon tiring effect was done by lOGS04j . predicting that while 
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the mass loss rate would remain the same, the velocity of the wind would drop. The new 
equation for the wind velocity becomes: 

^/^^) ^ ^cff[ro(x)] Y\l-m{w + x)]-l. (19) 

K 

Note that Fg is no longer a constant but decreases with the radius and the velocity. The 
photon tiring number can be found by combining equations fllSI) and fllSp . For the special 
case where a-p = 1/2 this implies that, 

m = 0.13 — — — , (20) 

r]oa2o Mq 

with the sound speed in units of 20 km/s. The implications of this equation are quite 
clear. The relative effect of photon tiring increases with Fg and decreases with the porosity 
length. This is to be expected, as a larger Fg means an increase in mass loss rate, and 
therefore an increase in the amount of energy necessary to lift the material. A larger rjQ 
means a decrease in coupling between radiation and matter, such that the effect of photon 
tiring diminishes. The same is true for an increase in Op, so we can expect the effect of 
photon tiring to diminish with higher a^. 



6.2 Numerical method 

We calculate the effect of photon tiring by calculating the work integral along the radial 
gridline and subtracting the result from the total luminosity of the star. This provides us 
with the luminosity available at each radial grid point, which can then be used to accelerate 
the wind during the next time step. 
The work integral is given by: 

m{r) fifraddr, 

-'\r (21) 
JR* 

Therefore, the effective luminosity at a given radius equals: 



pr 
JR* 



(22) 



(Ivan Marie et al.ll2008al ). Using this value for the luminosity, rather than the stellar lumi- 
nosity used in the simulations described in ||5], a new set of simulations was made for the 
same range of Fe values and ap = 0.5. All other parameters were kept the same as in ^ 
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Figure 5. The effect of photon tiring on the wind velocity (I). Clearly, the effect is much larger for larger Fg. 

6.3 Numerical results 

The effect of photon tiring is demonstrated in figure [5l For Fe = 2, the effect is small, 
whereas for Fg = 6 the effect is very large. Also, for the larger value of Fg, the velocity 
decreases beyond a certain radius. This corresponds to the radius where the luminosity falls 
below the Eddington luminosity. Note that the mass loss rate is independent of the tiring 
parameter and will therefore not change if photon tiring is included. This implies that the 
density of the wind increases, since the velocity decreases due to the photon tiring effect. It 
therefore also implies that the optical depth of the wind is higher as well. 

In figure El one can see a comparison between the values of the terminal velocity obtained 
with and without photon tiring. The influence of photon tiring clearly increases with Fg. For 
large Fg, the terminal velocity itself actually decreases. It is clear from these results that for 
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Figure 6. The effect of photon tiring on the wind velocity (II). This graph depicts the terminal velocity as a function of Fe 
for simulations with and without photon tiring included. The effect clearly increases with Fe. 



larger values of Fg the photon tiring effect must be included in the numerical simulations, 
or the results will quickly become unphysical. 



7 DISCUSSION 

A series of numerical simulations was carried out in order to test the analyt ical app roxi 



OGS04. The 



mations for the porosity length formalism of continuum driving, published by 
numerical results coincide well with the analytical results and demonstrate that this mecha- 
nism allows for powerful radiation-driven winds. This effect can explain the mass loss rates 
and wind velocities observed in Luminous Blue Variables such as rj Carinae. 

The simulations also confirm that the photon tiring effect plays an important role in 
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continuum driving as it places an upper limit on the mass loss rate. The effects of actually 
crossing the photon tiring limit have not yet been explored. This situation is much more 
complicated, since the simulations will no longer be able to reach a steady state solution. 
Ideally, such simulations should be done in two, or even three dimensions, to investigate 
the effect of interactions between different layers of the stellar wind as they move back and 
forth. Note that at this point we cannot predict how the star itself would react to such a 
situation. All our simulations have been done under the assumption that stellar parameters 
do not change significantly over time. It is possible that conditions in the outer layers of the 
star would change to reduce the driving force. 

Porosity reduced continuum driv ing caii also be important for the winds of other s uper- 



Eddington objects such as N ovae fjShaviv 



transients like M85OT2006 
In a companion paper, ( 



2001bl) . accretion disks (iBegelman 



Kulkarni et al. 



20061) and 



2007h. 



van Marie. Owocki fc Shaviv 



2008bl ). we explore the situation 



where the star exceeds the photon tiring limit. We also intend eventually to generalize the 
simulations to multiple dimensions and explore the influence of stellar rotation on continuum- 
driven winds. 
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